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Abstract 

This report describes some of the work done in a study of 
the flowfield produced by tangential leading-edge blowing on 
a 60-degree delta wing. The flow is investigated 
computationally by solving the Thin-Layer Navier-Stokes 
equations. Steady-state flowfields are calculated for 
various angles of attack and yaw, with and without the 
presence of tangential leading-edge blowing. The 
effectiveness of blowing as a rolling moment control 
mechanism to extend the envelope of controllability is 
illustrated at pre- and post-stall angles of attack. The 
numerical grid is generated using algebraic grid generation 
and various interpolation and blending techniques. The jet 
emanates from a slot with linearly-varying thickness and is 
introduced into the flowfield using the concept of an 
actuator plane, thereby not requiring resolution of the jet 
slot geometry. The Baldwin-Lomax algebraic turbulence model 
is used to provide turbulent closure. The computational 
results are compared with those of experiments. 
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Nomenclature 


C L wing lift coefficient 

C N wing normal force coefficient 

c Mx wing rolling moment coefficient 
Cjijy wing pitching moment coefficient 
blowing momentum coefficient 

a angle of attack 

J3 yaw angle 

freestream Mach number 
S ref wing reference area 

Re Reynolds number based on wing chord 

Uj jet velocity 

mj jet mass flow rate 

dynamic pressure 
AR aspect ratio 

Cp Jt „ stagnation pressure coefficient 



1. Introduction 


The extension of the high angle of attack (a) regime of 
delta wing aircraft has been of considerable interest in 
.recent years. It is well known that the flow at high a over 
delta wings is dominated by two leading-edge vortices. 
Previous studies have shown the occurrence of vortex 
asymmetry and burst, the first of which leads to the wing 
rock phenomenon. Both vortex asymmetry and burst cause non- 
linear and often unpredictable aerodynamic forces and 
moments. Since conventional control surfaces are ineffective 
at high a, it is desirable to have a mechanism which would 
control the strength and location of these leading-edge 
vortices. 

By modifying the vortical flowfield, tangential leading-edge 
blowing (TLEB) has been shown experimentally and 
computationally to be a viable option in providing roll 
control for delta wings at high a as in Refs. 1 and 2. These 
studies were for a half-span wing assuming symmetry. More 
recently, experimental studies were performed using a full- 
span wing. ' The effectiveness of TLEB was again 
demonstrated but there is coupling between the blown and 
unblown sides for pre-stall angles of attack. Data from 
these studies include force and moment coefficients as well 
as surface pressures and laser sheet flow visualization at 
one downstream station. 

Previous computational studies of the subsonic flow over 
delta wings at high alpha have mainly concentrated on wings 
with sharp leading edges 10-16 where separation is generally 
fixed. A round leading edge necessitates the use of a 
viscous method in order to determine the exact location of 
the primary separation line. The accuracy of the location of 
the separation line plays an important role in the 
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trajectory of the leading-edge vortex, and therefore on the 
aerodynamics of the flow. 

The current report describes part of a computational study 
solving the thin-layer form of the Navier-Stokes equations 
for the flow about the same full-span model as used in the 
aforementioned experiments. 3 ' 4 The numerical grid is 
generated algebraically using various blending and 
interpolation functions. The jet slot is not modelled 
explicitly in the grid (as by Yeh et al 2 for a half-span 
delta wing) , but is introduced using the actuator plane 
concept as used by Tavella et al 5 and Font et al. 6 The flow 
is solved using the F3D code, which has been successfully 
applied to high a flow. 5,6 ' 7 The Baldwin-Lomax half-equation 
algebraic turbulence model, with modifications suggested by 

9 

Degani and Schiff, is used to provide turbulent closure. 

The objectives of this study are: 

(i) to validate the current approach by comparing 
computational results with that of experiment, 

(ii) to provide a detailed description of the flowfield and 
aerodynamics governing the flow at angle of attack with 
and without the presence of blowing, and 

(iii) to show the effectiveness of TLEB as a roll control 
mechanism at high angles of attack. 

2. Geometry 

The delta wing used for this numerical experiment has a 
leading-edge sweep angle of 60° and is given in Fig. 1. The 
geometry is the same as that used for the experiments 

rj i 

reported by Celik et al and Wood et al. The leading edges 
are round and the top and bottom surfaces of the wing are 
parallel with a constant thickness (6% at the root chord) . 
The nose of the wing is generated by continuing the leading 
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edges to the apex. The result of this is a sharp nose when 
viewed from above but an elliptically shaped nose when 
viewed from the side. The wing is cropped to form a constant 
span tail or flap which is angled on the lower surface. The 
jet slots for blowing are situated on the wing leading edges 
from the nose-body junction to just short of the beginning 
of the wing tail. 

3. Grid Generation 

Due to the relatively simple grid topology required for this 
geometry, the grid about the delta wing is generated 
algebraically. Algebraic grid generation is efficient and 
therefore lends itself to grid refinement studies. The grid 
is generated in 2-D sections and then stacked together. 
Several post-processing steps are then performed, mainly 
through interpolation, to redistribute points in a desirable 
3-D fashion. Vinokur stretching is used in the 
redistribution, 

3.1 Topology 

The grid topology is C-O. There are three sections to the 
original grid. The first is the nose which has a spherical 
cap. The generation of this cap is done in 2-D semi-circular 
disks which are stacked in the circumferential direction. 
The second section is the body of the wing where 2-D 
downstream cross-sections or cross-flow planes are 
generated. In cross-section the leading edge is elliptical 
since it is circular normal to the leading edge. The last 
section is the tail-wake region where 2-D downstream 
sections are again generated. Cross-sections of the tail are 
rectangular due to the cropped wing. A zero-thickness wake 
with flow-through conditions is added to the wing to allow 
for closure of the computational domain. 
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Stretching is employed in the downstream direction on the 
wing body and in the wake. The grid extends six chord 
lengths in front and nine chord lengths behind the wing. 

3.2 Orthogonality 

1 8 

As suggested by Tavella, an algebraic grid can be improved 
by post-processing, particularly by blending onto desired 
functions close to a body surface. This is implemented by 
projecting the arc lengths of the existing grid at a surface 
onto a normal to the surface, and then blending the existing 
grid and the normal by using an inverse tangent function 
combination. The number of intervals blended, the grid 
location where the two functions being blended coincide 
(transition index) , and the power of the blending functions 
need to be specified. The blending function is given by the 
following equation : 

x(Q = x,(0.f, + x 2 (0.f 2 


where 



Lt = transition index 
p = power 

x = generic coordinate 
( = index 


The surface normal is determined by defining two body- 
tangent vectors at each surface grid point. The cross 
product of these two vectors divided by its length gives the 
unit normal, n. From the unit normal, we can calculate a 
projection matrix, P=n T n/n*n, which when post-multiplied by 
any vector would give the projection of that vector along 
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the unit normal. A cross-section of the grid shown in Fig. 2 
illustrates the result of blending. 

3.3 Zonal arrangement 

The implementation of the flow solver requires multiple 
zones of a certain maximum size because of memory 
limitations. To accommodate this, the grid is divided into 
fourteen zones, half of which are shown in Fig. 3. No 
interpolation is required at the zonal interfaces. Symmetric 
zones are used to save on secondary storage because the 
metrics are identical for both sides of the symmetry plane. 6 
Note also that zones are separated circumferentially at the 
blowing slots for the implementation of the actuator plane 
as an inter-zonal boundary condition. The grid is 80x137x50 
or 548000 points. For the generation of a grid of this size, 
the total Cray Y-MP CPU time required is 23 sec. 

4. Numerical method 

The flow is solved for using a version of the F3D code. F3D 

contains the algorithm developed by Steger et al. ' The 

implicit two-factor algorithm uses flux-vector splitting to 

upwind difference the convection terms in the streamwise 

direction, while retaining central differencing in the other 

directions. This scheme is unconditionally stable for the 

model wave equation and can have natural numerical 

dissipation and better stability properties than fully 

1 9 

central difference algorithms. 

Due to the high Reynolds number flow in this study, the 

g 

Baldwin-Lomax turbulence model, with modifications for 
cross-flow separation suggested by Degani and Schiff, 9 is 
used. Values for y + are consistently less than 4 throughout 
the grid. The wake is assumed to be laminar. The flow on the 
wing is treated as fully turbulent. 
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The convergence rate is relatively slow because of the low 
Mach number used (0.2), the rounded leading-edge wing (which 
means that the primary separation line takes time to be 
established in a steady-state calculation) , and due to high 
values of the grid transformation Jacobian in the stretched 
nose grid. 

5. Actuator plane concept 

The difference between using the actual geometry and using 
the actuator plane concept is shown in Fig. 4. Note that the 
actuator plane is, in effect, a discontinuity imposed as a 
boundary condition in the flow at an inter-zonal boundary. 
This concept has been used before on an ogive cylinder 6 and 
the F-18 forebody. 5 In the latter study, a comparison was 
made between resolving and not resolving the jet slot. The 
differences in results were found to be minimal for a delta 
wing test case. 

The jet Mach number is calculated assuming incompressibility 
or jet exit pressure equal to that of freestream. The 
boundary condition at the jet interface of the zone below 
the jet is obtained by extrapolation from the interior 
points . 

6. Results and Discussion 

All cases were run at M^ = 0.2 and Re=1.2xl0 6 based on the 
wing chord. 
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6.1 No blowing cases 

6 . 1.1 a=30°, ( 3 = 0 ° 

Figures 5, 6 and 7 show the stagnation pressure 

contours, simulated surface flow patterns and particle 
traces respectively for this case. Note that at this a 
the vortex stays unburst until it reaches the tail 
portion of the wing. This is evident from the sudden 
thickening of the core in both the particle traces and 
in the Cp stag contours. The primary and secondary 
vortices can also be seen on the Cp stag contours (Fig. 
5) . From Fig. 6, the primary separation line extends 
along the leading edge, and moves outboard. This 
variation in separation location is the main difference 
between this wing and sharp leading-edge wings where 
the primary separation line is fixed at the leading 
edge itself. The location of the primary separation 
line at the nose can be connected to a delay in the 
formation of the secondary line. An interesting feature 
to note is the formation of another vortex with 
rotation in the same sense as that of the primary one 
further downstream. Its development is then altered by 
the forced separation at the rectangular cross-section 
tail. Whether a whole series of vortices will be 
produced on an infinitely long wing is an intriguing 
possibility. A determining factor could be whether the 
primary leading-edge vortex feeding sheet stays 
attached to the primary separation line further 
downstream. 

Figure 8 gives the calculated spanwise pressure 
coefficient distributions at five downstream stations. 
Here again the formation of the vortex is evident. A 
signature of a rounded leading-edge wing is clearly 
visible in the form of the wide leading-edge suction 
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peak. For sharp leading-edge wings this peak is often 
only a spike, thereby not contributing significantly to 
the lift. 

6.1.2 a=50°, p=0° 


Figures 9, 10 and 11 show the stagnation pressure 

contours, simulated surface flow and particle traces 
for this case. At this a vortex breakdown is observed 
at 30% root-chord. The breakdown can be described as 
bubble type rather than spiral because of the fixed 
envelope which the streamlines and contours form and 
the extension of the secondary separation line 
downstream. Figure 10 shows that the primary separation 
line weakens further downstream because the flow at 
this a has a dominant cross-flow component. The 
surface flow patterns on the tail show a saddle point 
under the stagnated flow of the burst vortex. The 
vortex trajectory remains close to the wing due to the 
suction on the vortex and then follows the freestream 
direction in the wake. 

Figure 12 gives the calculated pressure coefficient 
distributions at the same five downstream stations for 
this case. The footprint of a burst vortex in this case 
is a flattening of the suction pressure peak in the 
downstream locations. The leading-edge suction peak is 
also thinner than in the a=30° case because of the 
more bluff body-type flow. 

6.1.3<x=30°, P=10° 

To test the code and the blowing mechanism for 
initially asymmetric flows, and to get at least a crude 
approximation of a real-life maneuver, a number of 
cases are calculated for a wing in side-slip. 
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Figure 13 shows the particle traces from a preliminary 
computation. The most notable feature of this yawed 
case is the intensification of the breakdown on the 
windward side and the absence of breakdown on the 
leeward side. From the pressure coefficient 
distribution in Fig. 14 it can be seen that the 
windward side vortex is stronger as expected. 

6.1.4 a=50°, (3=10° 


Figure 15 gives a preliminary plot of the pressure 
distributions for this case. The same features of 
asymmetry noted above are again evident. 

6.2 Blowing cases 


The blowing momentum coefficient is defined by 


C 

n 


1A 


<U5 


ref 


6.2.1 a = 30°, |3 = 0° 

A range of blowing coefficients was used to effectively 

show the influence of blowing on the aerodynamics. From 
the experiments, a value of of about 0.04 marks the 

boundary of the useful blowing region. Increasing the 
blowing intensity further produced no extra rolling 
moment at this pre-stall angle of attack. In order to 

try and capture this non-linear behaviour, a range of 
from 0 to 0.06 was chosen in increments of 0.01, 

The effect of tangential leading-edge blowing on the 
flow at this (X for three blowing coefficients is shown 
in the particle traces of Figs. 16, 17 and 18. The 
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blowing slot is indicated on the pilot's right with a 

thick line on the side of the wing. Note how an 

increase in blowing intensity effectively moves the 
leading-edge vortex inboard, until at some blowing 
coefficient the flow becomes completely attached. Note 
also, more importantly, that at this a, the flow on 

the unblown side remains essentially unchanged by the 
blowing, and is therefore uncoupled. At lower 

values, the leading-edge vortex on the blown side is 
not as focused or strong as on the unblown one, and 

breakdown seems to occur sooner. 

Figures 19, 20 and 21 show the stagnation pressure 

contours for three values of C^. The contours show how 

the vortex is flattened compared to the unblown side 
until it is eventually 'blown away' . 

Figures 22, 23 and 24 show the simulated surface flow 

contours. Here the effectiveness of blowing in moving 
the primary separation line is evident. Note that at 
the higher blowing coefficients the surface flow goes 
over the symmetry line, but does not significantly 
affect the unblown side flow. 

The calculated spanwise pressure coefficients at one 
downstream station as a function of increased blowing 
intensity is shown in Fig. 25. The uncoupling is 
evident in the collapse of the curves on the unblown 
side. On the blown side, the re-attachment of the 
leading-edge vortex with increased blowing can be seen. 
The lift contribution of the leading-edge suction peak 
also increases with increased blowing. This will later 
be seen to contribute to a blown-side-up rolling 


moment . 
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6.3 Comparison with experiment 

Before comparisons are made, it should be noted that the 
experiments were run at M^O.058, Re=4xl0 5 and 0.12, 

ke=8xl0 5 for all a and P cases, and that significant wind 

4 

tunnel wall interference at high a is noted due to a 
relatively small test section. It should be pointed out, 
however, that the aim of the experiment was to show the 
effect of blowing, and therefore the difference in 
coefficient values between the blown and unblown cases, 
rather than the absolute values of the coefficients. 

4 

Reservations are also expressed as to the accuracy of the 
stall angle and burst locations in view of comparison with 
other similar experiments. To summarize, comparison studies 
between computation and experiment should only be made if 
the difference between blown and unblown coefficients are 
used. Before comparing specific quantities, it should be 
mentioned that an angle of attack correction is expected to 
be justified. Specifically, it is suspected that the 
effective angle of attack in the experiment is higher than 
the geometric angle of attack due to wall interference 
effects. This is being investigated carefully. 

6.3.1 Comparison of pressure coefficients 

The difference in pressure coefficients when blowing is 
applied is used as a measure of comparing with 
experiment. ACp is defined as -Cp_unblown - (- 

Cp_blown) . This means that a positive ACp would 
correspond to a reduction in suction peak. Figures 26, 
27 and 28 show how this measures changes with blowing 
coefficient for a = 30°. The reduction in the pressure 

peak agrees reasonably well with experiment. The only 
real discrepancy is at low where the experimental 

difference goes slightly negative. This implies that in 
the experiment the leading-edge vortex spread out and 
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moved inboard as its strength was reduced by the 
blowing . 

6.3.2 Comparison of rolling moment coefficient 

Figure 29 gives the rolling moment comparison for the 
a = 30° case. Note how the saturation effect is 

captured as well as the saturation level. This confirms 
the notion m the experimental study that blowing for 
a pre-stall angle of attack is only useful in a limited 
range. Since the rolling moment is zero at zero 
blowing, the figure actually shows the difference in 
rolling moment coefficient between blowing and no 
blowing . 

6.3.3 Comparison of normal force coefficient 

The difference in normal force coefficients as blowing 
is applied is shown for the a = 30° case in Fig. 30. 
The difference is calculated the same way as above for 
the pressure coefficient. The same trend is observed in 
the computations as in the experiment, but there is a 
distinct difference in the turning point where after a 
reduction in normal force with increased blowing, the 
normal force increases again. One explanation for this 
discrepancy could be that in the computations, the 
calculated flow is further from the stall point as in 

4 

the experiment , where it was determined to be at a = 
40°. In other words, an angle of attack correction 
should be made as was mentioned earlier. The reason for 
this can be seen on the same plot where the a = 20° 

normal force coefficient differences are shown. The 
turning point is at a lower value for this angle of 

attack further away from the stall point. Another 
supporting fact is that the absolute normal force 
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coefficients are much higher in the experiments at 30° 
than in the computations. More computations are being 
planned to determine the stall angle for the 
computational configuration. Unfortunately the 
computations are too expensive to allow for a whole 
range of angle of attacks needed to determine the 
computational stall point. 

Finally, a comparison is made for the absolute normal 
force coefficient values. Figure 31 shows the 

computational C L values for the a=30°,p=0° and 

<x=50°,P=0° cases as compared to the Leading-Edge 

• 21 
Suction Analogy theory of Polhamus for incompressible 

flow and experiments 16 ' 22,23 ' 24 for sharp leading-edge 

wings at comparable AR. The tail on the wing under 

study gives it an effective AR=1.224. The Polhamus 

analogy is for sharp leading-edge and pure delta wings 

and is therefore only useful in quantitative 

comparisons. The theory is also only expected to be 

accurate for a's before breakdown occurs. 21 

Unfortunately, only normal force coefficients and no 

lift or drag coefficients were measured in the 

. 3 4 

experiments, ' so a direct comparison of lift is not 
possible . 

22 

Soltani et al showed that an increase in Reynolds 

number decreases the vortex lift at a given a. This 
could partly explain why the normal force coefficient 
values of the computations are below that of 
experiment . 

6.3.4 Comparison of pitching moment coefficient 

Here, we see essentially the same trend as for the 
normal force coefficient difference, as is shown in 
Fig. 32. The same arguments also hold. 
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7. Future work 

Grid resolution is another key factor in accurately 
capturing leading-edge vortices. In order to determine 
whether the lack of good agreement with experiment is due to 
experimental conditions (e.g. wind tunnel wall interference 
or lower freestream velocity) or due to numerical inaccuracy 
due to grid coarseness, a finer grid has been generated. The 
points are more concentrated in the vortex trajectory area. 
Computations with this configuration are being planned. 

Recent experimental results with a similar wing at the 
University of Bath has shed some light on the poor 
comparison of absolute values obtained above. Fig. 33 shows a 
comparison of the rolling moment coefficient contours in a- 
space for the Stanford and Bath experiments. The planform 

of the wings are identical, but the thicknesses are 6% and 
3.3% respectively. The blowing slots are different also, but 
most importantly, wind tunnel wall inteference is much less 
for the Bath experiments with a wing/tunnel area ratio of 
only 5% at a = 90°. The contours show the same basic 

control reversal, but for the Bath case it is confined to 
the stall angles of attack. For instance at a = 50° control 
reversal occured in the Stanford experiments, but not in the 
Bath experiments. Preliminary computational results at 50° 
also show no control reversal. Another point from Ref. 25 is 

A 

the values of C N obtained. In the Stanford experiments , the 
value at a = 30° was about 1.8, whereas for the Bath 

experiments it was about 1.4, another indication of the 
influence of wind tunnel walls. In the computations, the 
value is about 1.2. Computations at a = 40° are now being 
run to to try to determine the stall angle and whether 
control reversal is obtained. 
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Comparison with other Navier-Stokes codes such as TNS 26 29 
and different turbulence models are also being planned. 

8. Conclusions 

This report shows that a successful calculation of the 
complicated flowfield around a rounded leading-edge delta 
wing at angles of attack and yaw with and without tangential 
leading-edge blowing impossible with current resources. The 
results obtained thus far compare reasonably well with 

experiments if the difference in coefficients between 
unblown and blown cases is used. The effectiveness of 

blowing as a means to provide roll control at post-stall 

angles of attack is illustrated. Future work will extend the 
envelope of cases in order to fully validate the 
computational results with experiments. 
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Fig.l Perspective view of wing with blowing slot 



Fig. 2 Close-up of nose grid showing effect of blending 













Fig. 6 Computed surface flow patterns (a=30°,P=0°) 



















Fig. 13 Particle traces 




Fig. 14 
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Fig. 16 Particle traces (a=30° , (3=0° , Cn=0 . 01) 













Fig. 26 Comparison of change in pressure coefficient due to 
blowing between computation and experiment 
(a=30 ° , P=0 ° , Cn=0 .01) 










Fig. 28 Comparison of change in pressure coefficient due to 
blowing between computation and experiment 
(a=30°,P=0°,C^=0.03) 
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Fig. 29 Comparison of rolling moment coefficient as a function of 
blowing moment coefficient between computations and 
experiment (a=30°,B=0°) 
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Fig. 30 Comparison of change in normal force coefficient between 
computations and experiment (oc=30°,P=0°) 
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Fig. 31 Comparison of (a=30°,P=0°) and (a=50°,P=0°) lift 
coefficient with Polhamus theory 21 and experiment 
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Fig. 32 Comparison of change in pitching moment coefficient 
between computations and experiment (a=30°,P=0°) 






Fig. 33 Comparison of rolling moment coefficient contours in 

3 4 25 

a-C u space between Stanford ' and Bath experiments 
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